The classical mechanics of non-conservative systems 



Chad R. GallejQ 

Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, 91109 and 
Theoretical Astrophysics, California Institute of Technology, Pasadena, CA, 91125 

Hamilton's principle of stationary action lies at the foundation of theoretical physics and is applied 
in many other disciplines from pure mathematics to economics. Despite its utility, Hamilton's 
principle has a subtle pitfall that often goes unnoticed in physics: it is formulated as a boundary 
value problem in time but is used to derive equations of motion that are solved with initial data. 
This subtlety can have undesirable effects. I present a formulation of Hamilton's principle that is 
compatible with initial value problems. Remarkably, this leads to a natural formulation for the 
Lagrangian and Hamiltonian dynamics of generic non-conservative systems, thereby filling a long- 
standing gap in classical mechanics. Thus dissipative effects, for example, can be studied with new 
tools that may have application in a variety of disciplines. The new formalism is demonstrated by 
two examples of non-conservative systems: an object moving in a fluid with viscous drag forces and 
a harmonic oscillator coupled to a dissipative environment. 



Hamilton's principle of stationary action [T] is a corner- 
stone of physics and is the primary, formulaic way to de- 
rive equations of motion for many systems of varying de- 
grees of complexity - from the simple harmonic oscillator 
to supersymmetric gauge quantum field theories. Hamil- 
ton's principle relies on a Lagrangian or Hamiltonian for- 
mulation of a system, which account for conservative dy- 
namics but cannot describe generic non-conservative in- 
teractions. For simple dissipation forces local in time 
and linear in the velocities, one may use Rayleigh's dis- 
sipation function pQ. However, this function is not suf- 
ficiently comprehensive to describe systems with more 
general dissipative features like history-dependence, non- 
locality, and nonlinearity that can arise in open systems. 

The dynamical evolution and final configuration of 
non- conservative systems must be determined from ini- 
tial conditions. However, it seems under-appreciated that 
while initial data may be used to solve equations of mo- 
tion derived from Hamilton's principle, the latter is for- 
mulated with boundary conditions in time, not initial con- 
ditions. This observation may seem innocuous, and it 
usually is, except that this subtlety may manifest un- 
desirable features. Remarkably, resolving this subtlety 
opens the door to proper Lagrangian and Hamiltonian 
formulations of generic non- conservative systems. 

An illustrative example. To demonstrate the shortcoming 
of Hamilton's principle, consider a harmonic oscillator 
with amplitude q(t), mass to, and frequency oj coupled 
to an environment modeled by N harmonic oscillators 



with strength {A„}^ =1 , amplitudes {Qn(t)}n=i> masses 
{Mn}n=iJ an d frequencies {il n }^ =1 . The action for this 
system is 



S[q,{Qn}} = 
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Let us account for the effect of the {Q n } oscillators on 
q(t) by finding solutions only to the equations of motion 
for the {Q n } and inserting them back into ([!]), which is 
called integrating out. The resulting action, 
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is the effective action for q{t), though it is sometimes 
called a Fokker action [2]. The resulting dynamics for 
q(t) is open to exchange energy with the environment of 
the {Qn} oscillators and should thus be non-conservative. 
When TV is large this open subsystem ought to be dissi- 
pative. Q 



data) and G^(t 



(t) is a homogeneous solution (from initial 
— t') is the retarded Green function for 
the Q n oscillator. 

The last term in ^ involves two time integrals and the 
product q(t)q{t'). The latter is symmetric in t ■<-» t! and 
couples only to the time-symmetric part of the retarded 
Green function. Hence, the last term in ^ equals 
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when using the identity G^(t'-t) = G^ v (t-t'). Apply- 
ing Hamilton's principle to the effective action ([2| yields 
the equation of motion for q(t) 
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There are a couple of key points regarding Q. First, 
the second term on the right side depends on the ad- 
vanced Green function implying that solutions to Q 
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FIG. 1. Left: A cartoon of Hamilton's principle. Dashed 
lines denote the virtual displacements and the solid line the 
stationary path. Right: A cartoon of Hamilton's principle 
compatible with initial data (i.e., the final state is not fixed). 
In both cartoons, the arrows on the paths indicate the inte- 
gration direction for the line integral of the Lagrangian. 



do not evolve causally nor are specified by initial data 
alone. Second, the kernel of the integral in Q is sym- 
metric in time, which means that the integral describes 
conservative interactions between 9 and the {Q n } oscilla- 
tors. Consequently, Q does not account for dissipation, 
a time-asymmetric process, that should be present when 
N»l.' 

These undesirable features can be traced back to the 
very formulation of Hamilton's principle, which solves the 
problem: "Find the path q(t) passing through the given 
values qi at t = t% and 9/ at t = tf that makes the ac- 
tion stationary" (see left cartoon in Fig. [T]). Stated in 
this way, it is clear that Hamilton's principle is appropri- 
ate for systems satisfying boundary conditions in time, 
not initial conditions. According to Sturm-Liouville the- 
ory [3], the time- symmetric integration kernel in Q, 
which is a Green function itself, satisfies boundary con- 
ditions in time. Likewise, boundary conditions in time 
imply that the corresponding Green function is time- 
symmetric. Thus, there seems to be an intimate connec- 
tion in the variational calculus between boundary (ini- 
tial) conditions and conservative (non-conservative) dy- 
namics. 

In the remainder I formulate Hamilton's principle with 
initial conditions for general systems, report some conse- 
quences, and present some examples. 

Hamilton's principle with initial data. A hint for how to 
proceed comes from the previous example. The advanced 
Green function in ([3| and Q appears because the fac- 
tor q(t)q(t') couples only to the time- symmetric part of 
the retarded Green function. "Breaking" the symmetry 
by introducing two sets of variables, say q\ and q2, im- 
plies that q\(t)q 2 (t') will couple to the full retarded Green 
function, not just its time-symmetric part. Varying with 
respect to only q\ gives the correct force provided one 
sets 92 = qi after the variation [4]. This procedure is 
formalized for general systems below. 

Let q = {qi}fLi an d 9 = {<ji}iLi be a set of N gen- 
eralized coordinates and velocities of a general dynam- 



ical system. Formally, double both sets of quantities, 
q — > ((71,92) and 9 — > (91,92)- Parameterize both co- 
ordinate paths as qi, 2 (t,e) = qi. 2 (t,0) + erfi^f) where 
91,2(^7 0) are the coordinates of the two stationary paths, 
e <C 1, and rj\ y2 (t) are arbitrary virtual displacements. 
To ensure that enough conditions are given for vary- 
ing the action we require that: 1) 7/1,2 (^i) = and 2) 
9i(i/,e) = Q2(tf,e) and 9i(t/,e) = 92 (tf, e) for all e (the 
equality condition). The equality condition does not fix 
either value at the final time since the values they equal 
are not specified. After all variations are performed, both 
paths are set equal to each other and identified with the 
physical one, q(t) (the physical limit). See the right car- 
toon in Fig. [T] 

The action functional of 91 and 9*2 is defined here as 
the total line integral of the Lagrangian along both paths 
plus the line integral of a function K (discussed below) 
that depends on both paths {9* a }a=i ancl cannot generally 
be written as the difference of two potentials, 

S[q a ] = / (fti(9i,9i) + / dtL(q 2 ,q 2 ) + dtK(q a ,q a ,t) 



■1 

dt[L(q\,q\) - L(q 2: q 2 ) + K(q a ,q a ,t)]. (5) 



This action defines a new Lagrangian 

A(q a ,q a ) = £(9i,9i) - £(92,92) +K(q a ,q a ,t). (6) 

If K could be written as the difference of two potentials, 
V{q\) — ^(92), then it could be absorbed into the differ- 
ence of the Lagrangians in ([5]), leaving K zero [5]. Thus, 
a non-zero K describes generalized forces that are not 
derivable from a potential (i.e., non-conservative forces) 
and couples the two paths with each other. 

It is convenient, but not necessary, to make a change 
of variables to 9+ = (91 +q 2 )/2 and 9I =9*1—92 because 
9*_ — > and 9+ — > q in the physical limit. The conjugate 
momenta in the "±" variables, regarded as functions of 
the "±" coordinates and velocities, are found to be p± = 
dA/dq^, and the paths are parameterized as q±(t, e) — 
q±(t, 0) + eff±(t). The new action ^ is stationary under 
these variations if = [dS{q±]/de} e =o for all ff±, or 
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f) + (t)-p-(t) + ff-(t)-p+(t) 



(7) 



where the subscript denotes evaluation at e = and 

V+-P-= E^i V+iP-u etc. 

The equality condition requires 9*1 (tf , e) = q 2 (tf, e) and 
qi(tf,e) = q 2 (tf,e) so that ff_(tf) = andp_(ty) = 0. 
With ff±(ti) = it follows that the boundary terms in ^ 
all vanish. Thus, |7| is satisfied for any fj±(t) provided 
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that the two variables q± (t) solve 



dt 



dA 
dq± 



(8) 



Of course, one could have used the qi j2 coordinates in- 
stead to find dpi^/dt = dk/dq\ y 2 with pi t2 — dA/dqi :2 
regarded as a function of q\ l2 and qi i2 . 

In the physical limit ("p.L"), only the dp+/dt = 
dA/dq- equation in (18]) survives, yielding 



dp(q, q, 

dt 



dA 



dq- 



dL 



p.l. 



where the conjugate momenta are 



p{q,q) 



OA 
dq- 



dL 
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dK 



dK 
dq- 



p.i. 



p.l. 



(9) 



(10) 



When K = the generalized forces are derived from 
potentials and one recovers the usual Euler-Lagrange 
equations. A non-zero K can be regarded as a "non- 
conservative potential." 

In the physical limit, only the Eulcr-Lagrange equa- 
tion for the "+" variable survives. Hence, expanding the 
action in powers of q_, the equations of motion in ([9| 
and ( 10 1 also follow from the variational principle 







SS[q±] 



Sq-(t) 



(11) 
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Only terms in the new action ([5| that are perturbatively 
linear in <f_ contribute to physical forces. 

A new Hamiltonian A is defined by Legendre trans- 
forming the new Lagrangian, 

^(3l,2,Pi,2) =pWi -P2-Q2 - A(q ly2 ,q ly2 ) (12) 
= H(q\, pi) - H(q 2 ,p 2 ) - K(q lj2 , $1,2, t) 

where q\ and q 2 are now functions of their respective 



coordinates and momenta. Writing (12 1 in the "±" vari- 
ables gives 



A{q±,p±) =p+-q- +P--q+ - A(q ± ,q±) 



Both ( 12 ) and (|13[) can be written as 



A(Va,Pa) =p a q a - A(q a ,q a ) 



(13) 



(14) 



where a "metric" 



Cab 



is introduced to raise and lower the 



indices labeling the doubled variables: (1,2) in (12) and 
(+,— ) in (13). For the former c a b — diag(l,— 1) and 
for the latter c a b = offdiag(l, 1) so that p a q a = c ab p a qb 
(repeated indices are summed) where c ab is the inverse 
of c ao . The appearance of a metric in (14), the hint of 



"covariance" in (12) and (13), and the use of doubled 



variables suggest additional structure for the symplectic 
manifold [B]. 



Define new Poisson brackets by 

rr , n df dg df dg 
{{f'9}} =^-^-^-^. 



(15) 



which can be shown to satisfy Jacobi's identity. Then, 
Hamilton's equations follow by extremizing the action 
giving 

OA ■ dA 

q a =g^ = {{q a ,A}}, p a = -— = {{p a , A}}. (16) 

Note the index positions since they are raised and lowered 



by the metric c a &. In the physical limit, (16) becomes 



Hamilton's equations for a non-conservative system, 
• dH 



q = 



p 



dp 
dH 
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= {q,H}-[{{q_,K}}] pl , 
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(17) 



dK 
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{p,H}-[{{p-,K}}] p] 



The total time derivative of the energy function [T] 

dL 



h(q,q) 



dq 



L, 



(18) 



follows from the usual manipulations 1 , which here give 



dh 



dL 

"dt 



d dK dK 



dt dq- dq^ 



(19) 



p.l 



The amount of energy entering or leaving the system is 
determined by K when dL/dt = and can be found 
directly from the new Lagrangian. 

Example: Viscous drag forces. This new formalism can 
be used like the standard theory. Consider the following 
new Lagrangian, given in the "±" variables, 



A(x± ,x±) = m X- ■ 



(20) 



where n = 1 (linear) or 2 (nonlinear). The first term 
is the difference of the two kinetic energies (= mx? /2 — 
m2?2/2), and the second term is K. The new Lagrangian 
( 20 1 is unique up to terms nonlinear in x_ and its time 



derivatives, which don't contribute to physical forces (see 
(111). Using (111, or Q and (10), gives the equations of 
motion in the physical limit, mx l — —ax l \x\ n . For 
71 = 1 the force is proportional to —x l and for n = 2 it is 
proportional to — x l \x |. The former is Stokes' law for the 
drag force on a spherical object moving slowly through a 
viscous fluid and the latter is a nonlinear drag force for 
motions with large Reynolds number [7]. The key point 
is that these (nonlinear) equations for dissipative motion 
are derived from a (new) Lagrangian. 

To show that the resulting solutions from initial data 
are consistent with the new Hamilton's principle, it is 
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sufficient to consider slow motions (n = 1) for which 
the equations of motion are linear. In the "±" variables 
the new Euler-Lagrange equations are mx\ = —ax l ±. 
The physical limit implies that x + is determined by the 
physical initial data, x+(tj) = x*i and x+(ti) = while 
x_ is specified by final data, X-{tf) = = x_(tf), 
according to the equality condition. Because X- does 
not survive the physical limit, prescribing (trivial) data 
for X- at the final time is of no physical consequence. 
The resulting solutions are x~{t) = and x + {t) = 
Xi + mvi/a [1 — e -a(*-*i)/ m ], The former automatically 
imposes the physical limit so that x + (t) is the physically 
correct solution. Also, the new action is stationary for 
these solutions, as can be shown by direct substitution 
into 0. 



With K given by the second term of (20) it follows 
that h = mx 2 /2 and dh/dt = — a\x\ n+1 , which is pre- 
cisely the energy lost per unit time by the object through 
frictional forces from viscous drag. 

Example: Harmonic oscillator coupled to an environ- 
ment. Return to the example of a harmonic oscillator 
q coupled to an environment modeled by N harmonic 
oscillators {Q n } to demonstrate that the new framework 
gives the correct physical description for the open dynam- 
ics of q. Assume initial conditions q(U) — qt, q(ti) = Vi, 
Qn(U) = Qni, and Q n (ti) = V ni . The total system is 
closed implying that K = and the usual action is given 
by ([I]). Doubling the degrees of freedom, the new action 
is constructed as in ([5} but with K = 0. The effective 
action for the open dynamics of the q oscillator alone is 
obtained by integrating out the {Q n ±} variables, which 
satisfy M n Q n ± -\-M n Vl n Q n ± — X n q±. Subject to the 
initial conditions and the equality condition at the final 
time, the solutions are 



Qn+(t) = Q^(t) + 



K 

M, 



dt'Gtl{t-t')q + {t') (21) 



n Jti 



l - [t) ~ M n J t " 



(22) 



where QX' (t) = Q ni cos fi n (t - U) + V ni / Q„ sin 0„ (t - 1, ) 
is the homogeneous solution. The "+" variable evolves 
forward in time and satisfies the initial conditions while 
the "— " variable evolves backward in time because of 
the equality condition at the final time. This is a general 
feature of the "±" variables. 

Substituting these solutions into the action yields the 
effective action, 



A' 



Ses[q±] = / dtl m(q+q- - u 2 q + qJ) + q_ ^ X n Q^ 



CO 



A 2 rt 



/ 'dt'q-(t)G^(t-t')q + (t')\. (23) 



The factor q_(t)q+(t') in the last term is not symmetric 
in t «-» t' and couples to the full retarded Green function 
as opposed to just its time-symmetric piece as in Q. 
Applying (111 to (23 1, gives the equation of motion in 



the physical limit, 



N A 2 f tf 

mq + mcA = £ -f- / dt' (t - t')q{t') 

n=l Mn J *i 
N 



(24) 



Now, the Green function appearing in ( 24 ) is the retarded 



one, G { ™l(t~t') = 9(t-t')/n n smtt n (t-t'), and solutions 
to (24 1 evolve causally from only initial data. 



From S c s = J dt A e g, the effective Lagrangian is 
A c ff(g±,g±) = m(q^q + -uj 2 q_q + ) + q_F(t) 

+ [ dt'q-(t)~f(t-t')q+(t'), (25) 



which is non-local and history-dependent. Here, F(t) = 
23n=i AnQn (i) acts like an external force and j(t—t') = 
S^Li A 2 /(M„r2 n ) sinf2„(t — if). The last two terms in 



( |25| constitute an effective non-conservative potential, 
K c g , that is non-local in time and history-dependent. 



From (19), the energy function evolves as 



dh 
eft 



qF(t) 



dt' f(t - t')q(t') 



(26) 



where h — m(q 2 +uj 2 q 2 )/2 is the energy of the oscillator 
from (18). To see a familiar dissipation, choose trivial 
initial data for the {Q n } so that F(t) = and take each 
M n to be a constant, M. The coupling strengths {A„} 
are arbitrary so let A„ = Af2„ for A constant. Then, 

7(t - t') = (X 2 /M)d/dt'J2n=l COSn n(t ~ f). If N is 

so large that q essentially couples to a continuum of os- 
cillators then the summation becomes integration over 
cos ft (t — t'), which is a Dirac delta distribution (local 
in time). With these considerations, the frequency ui 2 is 



renormalized to w 2 cn = cu 2 — S(0)X 2 /(mM) and (26) be- 



comes dh/dt = — 7o<? 2 (£) for 70 = A 2 /(2M), which is the 
power lost by a damped, simple harmonic oscillator. 

Concluding remarks. The key aspects of this non- 
conservative mechanics are the formal doubling of vari- 
ables and the K function describing non-conservative 
forces and interactions. While I have focused on discrete 
mechanical systems, mostly for demonstrative purposes, 
it is equally applicable to continuum systems like field 
theories (see [8] for a non-trivial application) and elastic 
media. 

An open system, which can exchange energy by in- 
teraction with some other set of variables, will have a 
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non-vanishing K . In the drag force example, the under- 
lying variables that cause frictional forces on the object 
are not given or modeled. Hence, K is prescribed so that 
the resulting drag force is the desired one. However, if all 
the degrees of freedom of a closed system are known or 
modeled then K — until a suitable subset of those vari- 
ables are integrated out or coarse-grained. This results 
in an open system for the remaining variables and gener- 
ates a non-zero K as in the coupled oscillators example 
(see discussion after (25)). Thus, a non-zero K may be 



related to a measure of entropy for statistical systems. 
The formalism developed here can be quantized by re- 



placing the new Poisson brackets in ( 15 1 by commutators, 
yielding a canonical quantization for non- conservative 
quantum systems (in this manner, one can quantize both 
examples given earlier). Also, for closed systems (K = 0) 
the new action |5]) is the classical limit of the so-called 
"in-in" quantum theory [9]. To see this, consider the 
"in-in" generating functional, 



Z = Vq{Dq 2 exp 



dt(L[qi,q\] ~ L[q 2 ,q 2 ] 



Then, in the stationary phase approximation with h — > 0, 
the classical action, —ihhiZ, is easily shown to equal 
S[q a ] + 0(h). 

This new formulation of non-conservative systems may 
be useful for any method or technique that normally uses, 
or could benefit from using, Lagrangians and Hamiltoni- 
ans. These might include: developing partition functions 
for non-conservative statistical systems (see also [TU]). 
studying the phase space structure of nonlinear dissipa- 
tive dynamical systems, developing symplectic numerical 
integrators for systems with physical dissipation, provid- 
ing a more rigorous foundation for kinetic theory, and 
modeling nonconservative quantum mechanical systems, 
among other things. In extra physical degrees of 
freedom are introduced in a Lagrangian to parameter- 
ize absorptive processes within the paradigm of effective 
field theory (see also [HI 02] for recent applications). 
That work, in combination with the results presented 
here, may provide a powerful tool for studies of dissipa- 
tive systems. 
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